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Abstract 

In this paper we characterise the indeterminate case by the eigenvalues of the Hankel matrices 
being bounded below by a strictly positive constant. An explicit lower bound is given in terms of the 
orthonormal polynomials and we find expressions for this lower bound in a number of indeterminate 
moment problems. 

1 Introduction 

Let a be a positive measure on R with infinite support and finite moments of all orders 

(1.1) s n = s n (a) = J x n da(x). 

Jr. 

With a we associate the infinite Hankel matrix = {Hjk}, 

(1.2) Hjk = sj + k- 

Let H.n be the (iV+1) x (jV+1) matrix whose entries are Hjk, < j, k < JV. Since "Hat is positive definite, 
then all its eigenvalues are positive. The large N asymptotics of the smallest eigenvalue, denoted as 



A at, of the Hankel matrix TLn has been studied in papers by Szego [11], Widom and Wilf [jl3f| , Chen 



and Lawrence 0. See also the monograph by Wilf [14]. All the cases considered by these authors are 
determinate moment problems, and it was shown in each case that A at —* 0, and asymptotic results 
were obtained about how fast Aat tends to zero. 

The smallest eigenvalue can be obtained from the classical Rayleigh quotient: 

{N N N 
y~) Sj+kVjVk : v j = !> £ R, < j < N 
j=0 k=0 k=0 
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It follows that A at is a decreasing function of N. 

The main result of this paper is Theorem 1.1, which we state next. 



Theorem 1.1 The moment problem associated with the moments (1.1) is determinate if and only if 
limjv^oo A at = 0. 

We shall compare this result with a theorem of Hamburger |8|, Satz XXXI], cf. [jl], p. 83] or |1C, p. 70] 
Let be the minimum of the Hankel form T~Lm on the hyper-plane vq = 1, i.e. 

{N N \ 
Y s j+k VjV k : Vq = l,Vj € R, < j < N > . 
j=0 k=0 ) 

and let /i' N be the corresponding minimum for the moment sequence s' n = s n+ 2,n > 0, i.e. 



f/ N = min < 



mm 



N N 

£ s J -+fc+2^^ : t/ = 1, vj- € R, < j < JV 

j=0 fc=0 
' iV+1 iV+1 

X] s j+k VjV k : vq = 0,v± = 1, £ R, < j < N + 1 
j=0 fe=0 



The theorem of Hamburger can be stated that the moment problem is determinate if and only if at 
least one of the limits liuiTv^oo UN, limjv->oo Mjv are zero - 

It is clear from (1.3), (1.4) that fj,^ > Ajv and similarly fj! N > Ajv+i. From these inequalities and 
Hamburger's theorem, we obtain the "only if" statement in Theorem 1.1. The "if" statement will be 
proved by finding a positive lower bound for the eigenvalues An, cf. Theorem 1.2 below. 

We think that Theorem 1.1 has the advantage over the theorem of Hamburger that it involves only 
the moment sequence (s n ) and not the shifted sequence (s n +2)- In section 2 we give another proof of 
the "only if" statement to make the proof of Theorem 1.1 independent of Hamburgers theorem. 

If 

N 

(1.5) vrAr(x) := y~]vjX j , 

3=0 

then a simple calculation shows that 

(1.6) ^2 s j+k VjV k = / tt 2 n (x) da(x), 

0<j, k<N ^ E 



2 



and 



N ,2tt 

(i.7) E w *=/ 



2 d0 

2tt" 



We could also study the reciprocal of A^v given by 



(li 



A 



max 



2t 



— : ir N , J n N [x)da{x) = 1 



Let {pk} denote the orthonormal polynomials with respect to a, normalised so that pk has positive 
leading coefficient. 

We recall that the moment problem is indeterminate, cf. if and only if there exists a non-real 
number zo such that 

oo 

(1.9) ]T|pa^o)| 2 <oo. 

k=0 

In the indeterminate case the series in (1.9) actually converges for all zq in C, uniformly on compact 
sets. In the determinate case the series in (1.9) diverges for all non-real zq and also for all real numbers 
except the at most countably many points, where a has a positive discrete mass. 

If we expand the polynomial (1.5) as a linear combination of the orthonormal system 

N 

Kn{x) = ^CjPjix), 



then 





3=0 


/■27T 








Jo 




we have defined 



2 dd 

2^ 



^2 CjC k j Pj(e td )p k (e t9 )i^ = ^2 fcjkCjCk, 

» <~^* T ^° n 0<j, k<N 



0<j, k<N 



2tt 



(1.10) Kj k = J p 3 (e w )p k (e-^. 



Thus 



(1.11) - — = max ■ 

An 



N 



^! K-jkCjCk '■ Cj, C j — 1 



0<j,k<N j=0 

Since the eigenvalues of the matrix {Kjk)o<j,k<N are positive, and their sum is its trace, then 

-27T N 



(1.12) T -<'£ l K kk = / 5> fc (e ie ) 
*n ^ Jo ^1 



2 d0 

27r' 



Thus in the case of indeterminacy, 

i r 2 * ^ i , 4f)s 2 de 

< 1J »> ^< o EM<=' 9 )j 

JV J u fc=o 
which shows that 

(1.14) lim Aiv > ( 1^ -tKk — 1 , 
where 

(oo 
£|p*(e") 
fc=0 

We recall that for z € C \ R the number p(z)/\z — ~z\ is the radius of the Weyl circle at z. 
The above argument establishes the following result: 

Theorem 1.2 In the indeterminate case the smallest eigenvalue Xn of the Hankel matrixH^ is bounded 
below by the harmonic mean of the function p along the unit circle. 

We shall conclude this paper with examples, where we have calculated or estimated the quantity 



(1.16) po = 




This will be done for the moment problems associated with the Stieltjes-Wigert polynomials, cf. p| , []i~2"| , 
the Al-Salam-Carlitz polynomials p|, the symmetrized version of polynomials of Berg-Valent (||) lead- 
ing to a Freud-like weight ||, and the q -1 -Hermite polynomials of Ismail and Masson |J. 

If we introduce the coefficients of the orthonormal polynomials as 

k 

(1.17) p k (x) = J2^J xj 

3=0 

then 




and therefore 

oo k 

(Lis) p =EE^i- 

k=0 j=0 
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Another possibility for calculating pq is to use the entire functions B, D from the Nevanlinna matrix 
since it is well known that pi p. 123] 

(i.i B) ± [n(z) ? = mm_mm. 

k=0 

It follows that 

oo 

(1.20) Yl \Pk(e i6 )\ 2 = lm{B(e ie )D(e- i6 )}/smB. 

k=0 



2 Indeterminate Moment Problems 

In this section we shall give a proof of Theorem 1.1 which is independent of Hamburgers result. We 
have already established that if limjv-+oo A at = 0, then the problem is determinate. We shall next prove 
that if Atv > 7 for all N, where 7 > 0, then the problem is indeterminate. Since 1/Ajv < 1/7 for all N, 
and 1/Ajv is the biggest eigenvalue of the positive definite matrix (K,jk)o<j,k<Ni we S e t 

(2.1) KjkCjC k < ~J2\ Cj \ 2 , 

0<j,k<N ' j=0 

for all vectors (co, ■ ■ ■ , cjy) € C N+1 . If we consider an arbitrary complex polynomial p of degree < N 
written as p{x) = Ylk=o c kPk{x), the inequality (2.1) can be formulated 

p2tt m 1 r 

(2.2) J |p(^)| 2 ^<-yip(x)| 2 «fa(a:). 

Let now zq be an arbitrary non-real number in the open unit disc. By the Cauchy integral formula 

2?r 7 e M - z 

and therefore 

(2-3) |p(. )| 2 < 1^ W W )\ 2 f- r f-. 

Combined with (2.2) we see that there is a constant K such that for all complex polynomials p 
(2.4) \p(z )\ 2 <K f \p(x)\ 2 da(x), 
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where K = 1/(7(1 - \z \ 2 )). 

This inequality implies indeterminacy in the following way. Applying it to the polynomial 

N 



k=0 

we get 

N 

(2.5) ^|p fc (z )| 2 <^, 

fc=0 

and since N is arbitrary, indeterminacy follows. 

Remark. We see that the infinite positive definite matrix /Coo = {K-j,k} is bounded on I 2 if and 
only if A at > 7 for all N for some 7 > 0. Furthermore /Coo is of trace class if and only if po < 00 • The 
result of Theorem 1.1 can be reformulated to say that boundedness implies trace class for this family 
of operators. 



3 Examples 

We shall follow the notation and terminology for (/-special functions as those in Gasper and Rahman 
0- 

Example 1. The Stieltjes-Wigert Polynomials. 

These polynomials are orthonormal with respect to the weight function 
k 

(3.1) u(x) = -^exp(-k 2 (logx) 2 ), x > 0, 

V7T 



where k > is a positive parameter, cf. Q, j^]. They are given by 

1 n f \ 

(3.2) Pn (v) = (-l)^ + hq-,q)n~ 2 J2(l) ^H^)*. 

k=o W ? 

where we have defined q = exp{— (2A; 2 ) -1 }. 
It follows by (1.18) that 

< 3 - 3 > «. - E^fE*' (2l+1, (l) 

n=0 W.9jn fc=0 Wg 
00 00 n , v 2 
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Putting n = k + j, the inner sum is 



^ g k+j {q;q)k+j _ g k , , k+1 , 



j 



and hence 

oo 2(k+-) 2 

(3.4) Po = J2 q t \ 2<M0,g fc+1 ;g;g,g)- 

fc=0 w,9jfc 

We can obtain another expression for po- We apply the transformation |7|, (III. 5)] 

(3.5) 201 (a, b; c; q, z) = ^/ C ' 3<fe(a, c/b, 0; c, cg/6z; q, q) 

[bz/c; q)oc 

to see that 

'n\ 2 1 



We then find 

-. 00 fc 



fc=0 j=0 



A formula more general than (3.6) is 

Euj n fn\ 2 _ 1 (lo; q)jUJ 2k ~ j 



(g; ?)n w g ( w ; g)°o ^ (g; g)j(g; g)Lj 

and is stated in Pj. This more general identity also follows from (3.5) and the simple observation 

(g~ fc ;g)j (wjh (u;q)k-j f , u - 
iv/qy . 



(g 1 k /u;q)j (uj;q) k {q;q) k ^ j 

We have numerically computed the smallest eigenvalue of the Hankel matrix of various dimensions 
with the Stieltjes-Wigert weight from which we extrapolate to determine the smallest eigenvalue s = 
limjv^ooAAT of the infinite Hankel matrix for different values of q. This is then compared with the 
numerically computed lower bound I = 1/po- For q = ^ we have s = 0.3605 ... , I = 0.3435 .... The 
percentage error 100(s — l)/s is plotted for various values of q and is shown in figure 1. 



Example 2. Al-Salam-Carlitz polynomials. 
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Figure 1: Percentage error plotted for various values of q. 



The Al-Salam-Carlitz polynomials were introduced in Q. We consider the indeterminate polynomi- 
als Vn (x; q), where < q < 1 and q < a < l/q, cf. @. For the corresponding orthonormal polynomials 
{ Pk } we have by §, (4.24)] 



(3.8) ^|p fc (e 

k=0 

Therefore 

(3.9) po 

where 

(3.10) I n 



i9\\2 ('/<"'''• '/'' '" : '/)-x. , ,fl -,H . . ,H -iH 



{aq,q,q;q) c 



h(e w ,e^,aq;qe^, q e^;q,q/a). 



2n 00 



i6\\2 



k=0 



J (l-q n e ie )(l- 



2tt (aq,q,q;q)oo ^ " (q;q)n ^aJ 



dQ 



q n e~ ie ) 2tt 
(z, l/z; g)oo dz 
| =1 (1 - q n z)(l - q n /z) 2mz' 



Recall the Jacobi triple product identity @, 

oo 

(3.11) j(z) : = (9,2,1/2^)00 = ^ c k zk , 

t 

with 



k=— oo 



(3.12) c k = (-1) 



q k{k+l)/2 + g k{k-l)/2 
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Note that Ck = c-\~. 

Using the partial fraction decomposition 

q n q -n 1 _ g 2n 



l-q n z l-q~ n z (1 - q n z){z - q n ) 

we find by the residue theorem and the Jacobi triple product identity (3.11) that for n > 1, I n is given 
by 

,2nw \ j _ „nr>„„ I j( Z ) ,_nl „-nr,„„ / j( z ) 



{l-q 2n )(q;q)ooIn = q n Res ( YZ~^ ,Z = °J ~ g ~" Res _ q -n z ' z = 



q 2_^q c -k-i-q i^q c- k -i 

k=0 k=0 

oo 

[q - q ) ck, 



k=l 



while for n = 0, Iq is 



(q;q)ooio = [ T , — 4^ — tttt— = ~ Res ( tt-K^ ^ = °) 

7ui=i (1 -z)(z- 1) 2ttz V I 1 - *) 2 / 



-J> + i) c _ fe _ 1 = J>i)V (fc+1)/2 . 



fc=0 fc=0 



The conclusion is 



(3-13) Jo = rVD-ljV^ 2 , 

{q;q)oc k=0 

'» - rata*: ^ '""')■ »* L 



;o ° fc=i 



The above formulas can be further simplified. Using the Jacobi triple product identity (3.11) we find 
for integer values of n 



^ (-1) = o, 



k=— oo 

hence 



(3.14) ^(-l)VV") = -^(-l)V n V fe ^, n = 0,±l,.... 

A:=0 fc=l 

This analysis implies 

oo 

(3.15) (q; q)^! - q 2n ) I n = 2 ^(-l) fc ?(*) 



^nfc _ ? -nfe 



fc=l 
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Thus we have established the representation for n > 1 

(3.16) 4.-4C- D-D*" 1 «(')^^ 
(«; 9)oo ^ sm(nr) 



g = e 



It is clear that Jo is the limiting case of I n as n — > 0. The representation (3.16) indicates that I n is a 
theta function evaluated at the special point nr, hence we do not expect to find a closed form expression 
for I n . 

Example 3. Freud-like weight. 

In m Berg-Valent found the Nevanlinna matrix in the case of the indeterminate moment problem 
corresponding to a birth and death process with quartic rates. Later Chen and Ismail, cf. ||, considered 
the corresponding symmetrized moment problem, found the Nevanlinna matrix and observed that there 
are solutions which behave as the Freud weight exp(— -\/[xj"). I n particular they found the entire functions 



(3.17) B(z) 
where 

(3.18) 8i{z) 

(3.19) K 
Note that 



-5 (^ovW2), D(z) = -S 2 (K ^z/2), 

7T 



(-1)" 

Un + lV." 

r(i/4)r(5/4) 



An+l 



I = 0,1,2,3, 



(3.20) 5q(z) 

(3.21) 5 2 {z) 



cosh(zv7) + cos(z\/7) , 
cosh(zv7) — cos(z\fi) . 

If oj := exp(i7r/4) = (1 + i)/\/2, then a simple calculation shows that 



1 

27 



(3.22) 



B(x)D(y) - D{x)B{y) 



-2i r 



7T 



cos(w Kq\/x/2) cos(u]KQ\Jy /2) — cos{uj Ko\/y/2) cos(luKo y x/2) 



If x = e ld , and y = e 16 , then we linearise the products of cosines and find that the right-hand side of 
(3.22) is 



— - {cos[i^o(^ 



iO/2 



uj 3 e ld/2 + ue- ie/2 )/V2] + cos[K (oj 3 e 



cos + cue 



we 



-i6/2 



)/v / 2] 



we 



ie/2 
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We now combine the first and third terms, then combine the second and fourth terms and apply the 
addition theorem for trigonometric functions. We then see that the above is 

— {sinhfi^o cos(0/2)] smh[K sin(0/2)] + sm[K cos(0/2)] sm[K sin(0/2)]} . 

7T 

Thus we have proved that 

B {e id ) D (e- id ) - B (e~ ld ) D (e id ) 



(3.23) 



gi0 g— id 



' {sinh[K cos(0/2)] sinh[K sin(0/2)] + sin[if cos(0/2)] sm[K sin(0/2)]} . 



7r sm ( 



Thus in the case under consideration, after some straightforward calculations and the evaluation of a 
beta integral, we obtain 



(3.24) «=/'£K(e<f^ = S E 



(if /2) 2m+2n 



(2m + l)(2n + 1) m! n! (m + n)! 

n=U m,ri>0,m+n even 

Example 4. g _1 -Hermite polynomials. 

Ismail and Masson ||] proved that for this moment problem the functions B and D are given by 

(3.25) B(smh£) = ^ , D(sinhf) = r— fe 5 ,?e 5 ;g)°o, 

{q,q;q)oo (<?;<?)oo 

H (5.32)], [||, (5.36)]; respectively. Ismail and Masson also showed that |], (6.25)] 

(3.26) ^(sinhO^sinhr/) - J B(sinhr ? )D(sinh^) 

oo 

11 [1 - 2e" V sinh^ - e~ 2ri q 2n \ [l + 2eV +1 sinh£ - e 2v q 2n+2 } . 



-e ' 



%;<z)oc n=0 

We rewrite the infinite product as 

00 oo 

1 J a n b n = a a n 6 n _i, 

n=0 n=l 

and with sinh^ = e tS and sinhr/ = e~* e we get the following representation 
B(e i6 )D(e- ld ) - B(e- ld )D(e ie ) 



(3.27) 



gi6 g— i0 



1 oo 

7 — r- II t 1 + A( f - 2 i 2n + 4 ^ 3n + « 4n - 8 ^ 2n cos (20)] 

_. oo 



it 



Writing the infinite product as a power series in cos 2 9 and using 

we evaluate the integral of (3.27) with respect to d6/2n as 
(3.28) n= (^k±(^) V (-W+"*»> 



i — I — / / ' ' L I ■ ■ ■ I I — I — i I ' ' li 1 

fc=0 l<«i<...<n fe 

The formula (3.27) can be transformed further by putting cos 2 ip = — cos and p 2 = q, because then 

oo oo 

n + i n f + 4 ^ n c ° s °] = n t 1 + p4n - 2p2n c ° s ( 2 ^)] 

n=l n=l 

can be expressed by means of the theta function f?i(p; tp)- We find 

(3.29) H [(l+<f) 2 + V cosfl] = ___ ^(-^("^^(cosV), 

n=l W,9joo n=0 

where 

sin(2n + 1)V> 

C/ 2 n(cos^) = : — 

smip 

is the Chebyshev polynomial of the second kind given by 

(3.30) U 2n (x) = g + J) (-l)*z 2 (^)(l - *y. 

Similarly putting cos 2 (p = cos 8 we find 

oo 1 oo 

(3.31) Yl [(l+<f) 2 -4<f cosfl] =___^(-l)" g mt/2n(co S ^). 

n=l W'^°°„=0 

If we let [/* be the polynomial of degree n such that U2 n (x) = f/*(x 2 ), we get 



(3.32) B{e»)D{e-») - B(e~ ie )D(e ie ) 



e lu — e 

oo 



_L_ (-l) m ^ 1 ) + (™ +1 )^(-cos^)^(cos0). 



^)°° n, m =0 
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For non-negative integers k, l,r we have 



(3.33) C(k,l,r) := — [ (1 + cos 6) k (l - cos 0) 1 cos r dB 

27T Jo 



' ^ (-l) r B{k + + \)2Fi{k + \, -r; k + I + 1; 2), 



vr v ' K 2' T v 2 

which gives 



(3.34) ±-JJu* n (- CO se)U* m {cose) d6 



= 1 1 1 ( 2T+1 1 ) <- i) " +,c( *' '.»+•»-*-«). 

fc=0 z=o v 7 v 7 

Putting these formulas together we get a 5-fold sum for po- 
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